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Accretion disks around black holes radiate a significant fraction of the rest mass of the accreting 
material in the form of thermal radiation from within a few gravitational radii of the black hole 
(r < 2QGM/(?). In addition, the accreting matter may also be illuminated by hard X-rays from the 
surrounding plasma which adds fluorescent transition lines to the emission. This radiation is emitted 
by matter moving along geodesics in the metric, therefore the strong Doppler and gravitational 
redshifts observed in the emission encode information about the strong gravitational field around 
the black hole. In this paper the possibility of using the X-ray emission as a strong field test of 
General Relativity is explored by calculating the spectra for both the transition line and thermal 
emission from a thin accretion disk in a series of parametrically deformed Kerr metrics. In addition 
the possibility of constraining a number of known black hole spacetimes in alternative theories of 
gravity is considered. 


I. INTRODUCTION 

General Relativity (GR) has been extensively tested 
in the weak field regime. The famous classical tests of 
GR began before the theory was even fully formulated 
with the explanation of the anomalous precession of the 
perihelion of mercury [T]. Only a few years later fur¬ 
ther confirmation came with the first measurements of 
the gravitational deflection of light during the 1919 so¬ 
lar eclipse [2]. More recently, but no less famously, the 
confirmation by [5] that the loss of energy from the bi¬ 
nary pulsar PSR1913-I-16 was in the manner predicted 
by the theory won the authors the 1993 Nobel prize in 
physics [3]. Interspersed among these landmark experi¬ 
ments GR has been the subject of a constant stream of 
increasingly stringent tests; to date no deviations from 
the predictions of GR have been detected. For a review 
of experimental tests of GR see [5]. 

All of the tests mentioned above concern phenomena 
in weak gravitational fields; it has proved much more 
difficult to devise similarly stringent tests of GR in the 
strong gravitational field. This is due in part to the im¬ 
mense difficulty associated with measuring the motion of 
bodies on the small gravitational length scale associated 
with black holes (a few kilometers for a solar mass black 
hole) at astrophysical distances. If such measurements 
could be made and it were possible to accurately track the 
motion of a test particle around an astrophysical black 
hole it would become possible to place constraints on the 
strong field spacetime metric. It is a prediction of GR, 
together with the no-hair theorem, that the spacetime 
metric around an astrophysical black holes is described 
by the famous Kerr solution [6]. Alternative theories of 
gravity predict the existence of different solutions (ex¬ 
amples of known solutions considered in this paper in- 


* cjm96@ast.cam.ac.uk 
^ jrg23@ast.cam.ac.uk 


elude the Kehagias-Sfetsos solution [7] to Hofava gravity 
BM and the slowly rotating solution uni to dynamical 
Chern-Simons gravity m)- If a deviation from the Kerr 
solution was detected it could indicate a failure of GR in 
the strong field regime. Alternatively, it could indicate 
a failure of the no-hair theorem and raise the possibility 
of exotic compact objects within GR. It should be noted 
that there are other ways in which alternative theories of 
gravity may differ from GR besides changes to the metric 
around a black hole; for example, scalar-tensor theories 
of gravity (of which Brans-Dicke theory [T^] is perhaps 
the best known example) generically predict the existence 
of addition gravitational wave polarisation states [5]. In 
this paper however we restrict our attention to testing 
the spacetime metric around a black hole. 

Tests of the metric around a black hole may be pos¬ 
sible in the near future with gravitational wave obser¬ 
vations, such as those from ground-based detectors like 
advanced LIGO and advanced VIRGO mui. For ex¬ 
ample, the possibility of using a network of ground-based 
gravitational wave detectors observing the coalescences 
of binary neutron stars to constrain the deviation of a 
post-Newtonian coefficient from the predicted GR value 
was considered in m- Further into the future, low- 
frequency space-based gravitational wave observatories 
such as eLISA [TB] will offer the possibility of probing 
the gravitational field around supermassive black holes. 
The possibility of using observations of the inspiral of a 
stellar mass compact object into a supermassive black 
hole to constrain the anomolous quadrupole moment of 
the massive black hole was considered by and later 
by [18] . These gravitational wave observations are ideally 
suited to such tests as they offer the ability to accurately 
track the orbital phase of the two bodies all the way up 
to and including the merger. However, as gravitational 
wave observation are not yet available it is interesting to 
ask if strong field tests are possible using current electro¬ 
magnetic observations. 

Accretion disks provide an ideal candidate for such 
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tests, because X-ray irradiation of matter in the inner 
regions of the disk (r < 20M) imprints characteristic 
features upon the X-ray spectra, in particular the fluo¬ 
rescent Ka Iron line (rest energy 6.38keV). As this emis¬ 
sion originates from so close to the Black Hole (BH) it 
is strongly distorted by gravitational and Doppler shift¬ 
ing producing the characteristic ‘two-horned’ profiles ob¬ 
served in nature m- Furthermore the light, once emit¬ 
ted, is strongly gravitationally lensed as it propagates 
through the spacetime, altering the observed spectra and 
imprinting upon it extra information about the strong 
field metric. Indeed the large width of these lines is evi¬ 
dence for the existence of highly spinning black holes ED]; 
for a recent review of black hole spin measurements us¬ 
ing X-ray emission see [Hj. In addition, viscous torques 
present in the disk dissipate energy causing the material 
to gradually spiral inwards, this energy is radiated locally 
in the form of thermal emission. This thermal radiation 
is subjected to the same strong-held gravitational effects 
as the line emission, and so may also be used as a probe 
of the metric. 


The approach of this paper is to consider Iron line 
and thermal emission in a large class of parametrically 
deformed Kerr black holes. The deformed black hole 
spacetimes are referred to as “bumpy black holes” and 
the individual deformation parameters as “bumps”. The 
spacetimes all have the property that if all the bump pa¬ 
rameters are simultaneously set to zero then the Kerr 
solution is recovered. The advantage of this approach, 
compared with considering the emission in a small num¬ 
ber of known alternative black hole solutions, is that one 
is able to consider a wide range of different deformations 
simultaneously and identify particular “bumps” which 
are especially easy or difficult to constrain. Furthermore, 
even if a particular alternative black hole is not contained 
within this class it may be hoped that the disk spectrum 
in this new metric would have a significant overlap with 
the spectrum of a black hole in the class; and therefore a 
deviation from GR would still be detectable. The family 
of bumpy black holes used for this purpose was the class 
of metrics constructed by [52]; these have the property 
that the perturbed spacetimes possess a fourth constant 
of motion (analogous to the Carter constant [53] in Kerr 
spacetime). In addition to this very general familiy of 
bumpy black holes the results are put into context by 
comparing them with the bounds it is possible to place on 
some known black hole solutions; in particular both the 
linear m and quadratic |24] in spin solutions to weakly 
coupled dynamical Chern-Simons gravity, and the slowly 
rotating Kehagias-Sfetsos black hole [25] were considered. 


In this paper we begin in Sec. |Tl| by describing the 
family of bumpy black holes that were considered in this 
paper, we also describe a couple of specific examples of 
known black hole solutions in alternative theories of grav¬ 
ity which are also considered here. In Sec. Ill the neces¬ 
sary theory for calculating both the thermal and Iron 
line emission from an accretion disk in a general sta¬ 
tionary, axisymmetric spacetime is described. Sec. jlVj 


describes the methods used for data analysis is this pa¬ 
per and the method used for estimating the bounds it 
will be possible to place on the different deformation, or 
“bump”, parameters. The results are presented in Sec. 


VI and finally a discussion and concluding remarks are 


given in Sec. |VII[ Throughout this paper natural units, 
where G = c = ks = 1, are used. 


II. BUMPY BLACK HOLE SPACETIMES 

A metric that represents a rotating black hole is re¬ 
quired to be stationary and axisymmetric (i.e. there ex¬ 
ists a timelike Killing vector d/dt and a spacelike Killing 
vector d/d(j)), invariant under simultaneous inversion of 
the (j) and t coordinates, and reflection symmetric about 
the equatorial plane. A sufficiently general metric which 
captures all these properties is given by [5^ 

ds^ = gtt dt^-|-5t0dtd0-|-5rr dr^-l-gee d0^-|-g00 d(()^ , (1) 

where the metric coefficients depend only on the radial 
and polar coordinates r and 9. In fact Eq. [^ retains 
a considerable degree of gauge freedom, however, it is 
sufficient for our present purpose. 

As our interest here is in testing the hypothesis that 
the metric is the Kerr solution, and it is known that the 
Kerr solution is an excellent description in the weak field, 
it is natural to expand the metric as 

9tit^ = j (2) 

where e 1 and h^i, —>■ 0 as r —>■ oo. The background 
Kerr solution is given by, 

„Kerr _ (^ 2Mr\ „Kerr _ —2M^ar „;„2 q 

9tt - - , 9t<p - —p— sm 9 , 

^Kerr _ „Kerr _ „2 „Kerr _ a (o\ 

9rr - A ’ 9ee - P ^ 9^^ - sin 9, (3) 

for a BH with a mass M and dimensionless spin a, in 

which 

p2 _ j,2 -1- q 2^2 j,Qg2 Q ^ A = (l — -I- , (4) 

and = (r^-I-A sin^ d . 

From here on dimensionless units with M = 1 are used. 

Geodesic motion with four-momentum in a met¬ 
ric with the symmetries outlined thus far would posses 
three constants of motion; energy E = —p^{d/dt)^, z- 
component of angular momentum L = p^^{d/d(f>)^, and 
particle rest mass rri^ = —p^Pfj.. However the Kerr so¬ 
lution also possess an additional constant of motion re¬ 
lated to the existence of a second rank Killing tensor, 
C = p^p'^^niy. In place of C the Carter constant is often 
defined as Q = C — {L — aE)^ [53]. The existence of the 
addition constant of motion ensures that geodesic motion 
in the spacetime is separable and tri-periodic (i.e. with 
a well defined frequency associated with motion in each 
of the r, 9 and (p coordinates). Although it is not clear 
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that it should be required of any alternative to the Kerr 
solution to possess these properties, they are sufficiently 
appealing that it is worth considering the possibility care¬ 
fully. 

The most general form of metric perturbation, 
that could be added to the Kerr solution such that the 
resultant metric also admitted a second rank Killing ten¬ 
sor and hence a Carter-like constant (at least to 0{e^) 
in the perturbation) was considered by |22j . The metric 
perturbation was also required to tend to zero in the limit 
r —>■ oo faster than r~^ so that all weak field tests are still 
satisfied. A similar solution was also found earlier by m- 
Working in Boyer-Linquist-like coordinates [32] found a 
series of differential relations that must be satisfied by 
the metric perturbation components. The only non zero 
components of the perturbation are {/itt,/irr, ^ 00 }; 
in particular hgg vanishes. In subsequent work |28] ex¬ 
pansions for the metric perturbation components in pow¬ 
ers of 1 /r were found. 



where the coefficients „ are functions of 9 only. The 
leading order non-zero coeffeicients are given in Eq. 
and all the coefficients up to 0{l/r^) are reproduced in 
Appendix |A 1[ 

htt ,2 = 71,2 + 274,2 - 2073,1 sin^ 9 , 

hrr ,2 = — 71,2 ) 

ht ^,2 = -M sin^ 9 [ 73,3 -|- a ( 71,2 + 74 , 2 ) + a^ 73 ,i] , 
/i 00 ,o= 2M^ 073,1 sin'^ 9 . ( 6 ) 

Note that the leading order correction enters at a lower 
order in the component. The components htt and 
hrr are dimensionless, whilst ht,p has units of length and 
has units of length squared; reflecting the dimensions 
of the components of g^,^. 

At low order in the expansion the metric perturba¬ 
tion is fully characterised by a small number of coef¬ 
ficients, 7 i,j. Adopting the notation of [25| the met¬ 
ric perturbation up to 0(l/r^) is given by the four 
constants B 2 = ( 71 , 2 , 737 , 73 , 3 , 74 , 2 ) (see Eq. [^; up to 
0(l/r^) it is given by the 7 constants B2 U S 3 , where 
^3 = ( 71 , 3 , 73 , 4 , 74 , 3 ); up to C>(l/r'‘) it is given by the 10 
constants S 2 US 3 US 4 , where S 4 = ( 71 , 4 , 73 , 5 , 74 , 4 ); up to 
O (1/r®) it is given by the 13 constants S 2 U S 3 U S 4 U S 5 , 
where S 5 = (71,5,74,5,73,6)- 

Throughout this paper various Sjv limits will be re¬ 
ferred to, these correspond to setting all the parameters 
7 i,j to zero except for those in the set Sat, which are 
treated as independent free parameters. This greatly re¬ 
duces the number of free parameters and allows for a 
systematic way of examining perturbations to the Kerr 
metric. Because the Kerr metric is known to be an excel¬ 
lent approximation at large radii it is natural to expect 
any deviation from the Kerr solution to show up, initially 
at least, at lowest order in 1/r. This is the justification 


for examining each of the Bn limits separately. In gen¬ 
eral, it would be possible for any combination of the 77 ^ 
to be be none-zero. In this paper results are presented 
for the B 2 , Bz, B^ and B^ metrics. 


A. Known black hole solutions 


In addition to the very general family of perturbed 
metrics described above, it is also useful to consider some 
examples of known black hole solutions in alternative the¬ 
ories of gravity. In this section a few such solutions are 
listed. One of these spacetimes (the linear in spin solu¬ 
tion to dynamical Chern-Simons gravity. Sec. IIA2) is a 
member of the family of solutions described in Sec. ^ it 
can be obtained by a particular choice of the constants 
li,3- 


1. The Kehagias Sfetsos metric 

A spherically symmetric black hole solution to Hofava 
gravity [3 H] was found by [7|. This was generalised to 
a slowly rotating solution by |2l|. Accretion disk sig¬ 
natures for this type of black hole have been considered 
previously by |29j . The metric is 

dsKS= -/Ks(4')d<^ -h - — rw + dfl^- d^dt , 

jKsir) r 

where fKsir) = 1 -f - ^1+ , 

and dfl^ = d9^ + sin^ 0d(/^ . (7) 

In the limit w —>■ cxd the slowly rotating limit of the Kerr 
metric is recovered. In order to avoid a naked singularity 
at the origin an extra constraint is needed, > d. 

For the remainder of this paper w is made dimensionless 
by multiplying by M^, and we choose to work with the 
small parameter Y = I/(a;M^) ^ 1. For the remainder 
of this paper the metric in Eq. [^will be refered to as the 
KS metric. 


2. The Chern Simons metric linear in spin 

Dynamical Chern Simons (CS) modified gravity | 11 ] 
is a parity violating theory of gravity constructed by 
adding a Pontryagin invariant term to the action. As the 
Schwarzschild solution is spherically symmetric and has 
even parity it remains a solution to the modified CS field 
equations; however, the Kerr solution does not have even 
parity and fails to satisfy these equations. No complete 
rotating black hole solution is known in the CS theory, 
however perturbative solutions in the spin and CS cou¬ 
pling constant have been found analytically. The rapidly 
rotating case was considered numerically by [30] , here we 
focus on the slowly rotating solutions. 
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The slowly rotating black hole solution (linear in spin) 
to dynamical Chern-Simons gravity was found by [10) : 

5 1 / 12 27 \ 

d4si = dsLrr + ^ + 1^ j 

Accretion disk signatures for this type of black hole have 
been considered previously by m- For the remainder of 
this paper the metric in Eq. will be refered to as the 
CSl metric. 


differential equation governing the radial dependence of 


the flux to be integrated; see Eq. 23 


The model described above is often referred to as the 
standard relativistic model of accretion disks, and has 
been extensively studied for the Kerr geometry. Here it 
is necessary to retain sufficient generality to perform the 
calculations for any of the bumpy black holes described 
in Sec. El Since the disk is thin and assumed to lie in 
the equatorial plane we may switch from using Boyer- 
Lidquist polar-like coordinates to cylindrical-like coordi¬ 
nates. The metric near the equatorial plane is given by 


3. The Chern Simons metric quadratic in spin 

A slowly rotating black hole solution (quadratic in 
spin) to dynamical Chern-Simons gravity was found by 
[24] . As far as the authors are aware disk emission in this 
metric has not been considered before. The metric is 

dscs2 = + 5 dx'^dx^ (9) 

and expressions for the components 6 are given in 

Appendix |A 2| Eor the remainder of this paper the metric 
in Eq. j^will be referred to as the CS2 metric. 

Unlike the KS metric, both the CSl and CS2 metrics 
recover the Schwarzschild solution in the limit a —>■ 0. 
This is associated with the fact that dynamical Chern Si¬ 
mons modified gravity is a parity violating theory of grav¬ 
ity and hence the spherically symmetric Schwarzschild 
solution in GR, which is parity even, remains a solution 
of the modified theory. 


III. DISK EMISSION IN BUMPY BLACK 
HOLES SPACETIMES 

The theory of thin accretion disks in the steady state 
was developed by [32| in the Newtonian case and gen¬ 
eralised by [33] to the general relativistic case. It is as¬ 
sumed that the disk lies in the equatorial plane with ma¬ 
terial moving in (approximately) circular geodesic orbits, 
with a small radial accretion velocity superposed on the 
circular motion. In addition, it is assumed that the disk 
is thin (i.e. at a radius r the thickness 2h satishes h <C r; 
the disk height may depend on radius, h = /i(r), but for 
notational simplicity this dependence is suppressed here) 
and is in a steady state with a constant accretion rate 
Mq. The disk extends from the innermost stable circu¬ 
lar orbit (ISCO) of the spacetime to some finite outer 
radius (strictly, to be consistent with the steady state 
assumption, the disk must extend indefinitely; here the 
more conservative approach of leaving the outer radius as 
a free parameter is taken). When the accreting material 
reaches the ISCO it plunges quickly into the BH with¬ 
out radiating any further or interacting viscously with 
the material at larger radii. The fact that the mate¬ 
rial plunges quickly after crossing the ISCO provides the 
physical motivation for imposing a zero torque boundary 
condition at the inner edge of the disk. This allows the 


ds^ = gttdt‘^+2gt^dtd4>+grrdr‘^+g^^d(j)‘^+gzzdz‘^ , (10) 


where the metric components depend only on r; correc¬ 
tions that depend on z enter at second order and are 
neglected because the disc is assumed to be thin. By a 
rescaling of the z coordinate the metric component gzz 
may be set to unity without any loss of generality. 

The metric in Eqs. and W does not depend on the 
timelike or azimuthal coordinates, so the corresponding 
covariant components of the four-momentum are con¬ 
served, pt = —E and = L. Using the metric to raise 
the indices on these momentum components gives the 
first two geodesic equations in first order form, where a 
tilde denotes an orbital quantity per unit mass of the test 
particle. 


dt 

dr 

d(j) 

dr 


Egrpfp -\- Lgt(p 
2 ’ 
9t4> 9tt94>4> 

Egt(j> + Lgtt 
9t4> ~ Sttgtfxp 


( 11 ) 


where r is the proper time along the worldline of the 
orbiting test particle. 

Because the material in the disk is orbiting in the equa¬ 
torial plane the vertical component of its four-velocity 
vanishes, dz/dr = 0. Therefore, the third and final 
geodesic equation may be conveniently obtained from the 
normalisation condition of the four velocity; 


dx>^ dx’' _ a u _ . 

9f^i^ dr dr — ^ ^ ? 


( 12 ) 





VeB{r) 

9Tr 


(13) 


where. 


Ve{r) — ^ 


9t^ 


gttgtjxj) 


(14) 


Erom Eq. 13 it can be seen that the conditions for sta¬ 
ble circular orbits are I4fi(r) = VeS,r{r) = 0, where a 
comma in the subscript denotes a partial derivative with 
respect to all subsequent indices. These two conditions 
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yield expressions for the specific energy and angular mo¬ 
mentum per unit mass of a particle on a circular, equa¬ 
torial geodesic; 


9tt + 

9t(j> + gtjitj)^ 

\/— {gtt + ‘^gtip^ + 9 (p(p^'^) 


(15) 

(16) 


Combining Eqs. |11| with Eqs. |15| and |16| gives an expres¬ 
sion for the coordinate angular velocity of the particle, 



~9t(l>,r + 



94>4>,r 


9tt,r9(f>4>,r 


(17) 


In the metrics of interest in this paper (i.e. the Kerr 
metric and small perturbations from it) there exists an 
ISCO. The radius of the ISCO may be found by solving 
14 ff,rr(?’) = 0 simultaneously with Veff(r) = I4ff_r(''’) = 0. 
This equation admits the following analytic solution in 
the Kerr case [34] : 

= 1 + (1 - a")'/" ((1 + a)i/3 + (1 - a)i/3) , 

^2= + , 

= 3 + ^2 - ((3 - Zi)(3 + + 2Z2)f'^ . (18) 

Eor the general bumpy black holes discussed in Sec.jTTjthe 
radius for the ISCO must be found numerically. However, 
if the spacetime is characterised by a small deformation 
from Kerr, then the radius of the ISCO may be expanded 
perturbatively in the bump parameter, e, 


Tisco = r 


Kerr 


dVetf,. 

r-{r,e) 


de 

r=r?"J‘,£=0 

dVeff,,. 

r(r.e) 


dr 



+ C>(e2) . (19) 


If the metric is characterised by more than one bump 
parameter, e^, as is the case with the Bn spacetimes de¬ 
scribed in Sec. |lT| t hen the following replacement should 
be made in Eq. |19[ 


,ei=0 

( 20 ) 

Erom Eq. [m and the fact that the orbit is circular 
and equatorial, the four-velocity of a small fluid element 
of the disk orbiting the black hole is given by 


dVef 




de 


be=0 




rr ) 


dei 


dt dr d(j) dz 
dr ’ dr ’ dr ’ dr 


^ = (m*, m’', u'^, = u* (1,0, H, 0) , 


,( 21 ) 

where u* may be found from the normalisation condition 
on the four velocity in Eq. [T^ 




(dtt + ‘^^gtip + ^"^gipip) 


( 22 ) 


Viscous forces in the disk cause the orbiting material to 
dissipate energy and gradually move inwards to smaller 
radii. An expression for the radial dependence of the flux 
of energy from the disk was derived in |33| ; for complete¬ 
ness this derivation is reproduced in Appendix The 
resulting expression for the radial dependence of the flux 
is, 


F(r) = 


— ^, r 




/ T 

{e - riL^ L^r dr 


(23) 

where Mq is the accretion rate and g is the metric deter¬ 
minant. Since the disk is in thermodynamic equilibrium, 
and not heating up or cooling down, this flux of energy 
is radiated away from the disk in the form of a thermal 
distribution of photons. 


A. Line emission 

In addition to the flux of thermal radiation a higher en¬ 
ergy power-law component of hard X-ray photons is also 
observed [3S]. For supermassive black holes with cooler 
disks this component can in fact dominate over the ther¬ 
mal component. The hard X-ray power-law component 
is generally believed to be caused by inverse compton 
scattering of the thermal photons radiated from the disk 
by the hot surrounding corona. A certain fraction of the 
hard X-ray photons are radiated back towards the disk 
surface where, upon incidence, they produce fluorescent 
transition lines which may also be observed as a third 
component of the spectrum. One frequently observed 
line is that due to the Iron (Fe) Ka transition, which in 
its rest frame has an energy of 6.38 keV. The combined 
effects of gravitational redshifting and Doppler boosting 
broaden this line into the characteristic shapes shown in 

Fig. 12 

The observed line profile can be computed as follows. 
Let dFo(Ao) be the infinitesimal element of flux at energy 
Eq observed at infinity due to an infinitesimal element of 
the disk. If Io{Eo) is the specific intensity (i.e. the energy 
per unit time, per unit area, per unit spectral energy (or 
frequency), per unit solid angle) in the observer’s rest 
frame then 


dFo{Eo) = Io{Eo)dE , (24) 

where the element of the disk subtends solid angle ele¬ 
ment dE. If the disk has a specific intensity Ie{Ee) in its 
own rest frame, given the invariance of I/E^ along the 
photon’s worldline |3S], it follows that 

dFoiEo)=g^hiEe)dE, (25) 

where g = 1/(1 -f z) = Eq/E^ is the red-shift factor. For 
line emission the specific intensity may be approximated 
as a delta-function, /e(Ee) = e(re, g,f,)S{Ef, — Emt), where 
Eint is the energy of the transition line in its rest frame. 
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Redshift z 


1 0.33 0 - 0.2 - 0.33 
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Redshift z 

1 0.33 0 - 0.2 - 0.33 
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Redshift z 

1 0.33 0 - 0.2 - 0.33 
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Redshift z 
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— a=0.0 
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— a=0.9 
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-- r„„,=30 
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FIG. 1. The dependence of the Kerr Iron line profile on various parameters. All spectra were normalised to the same peak flux. 
Unless otherwise stated in the hgure all parameters were set to their fiducial values; a = 0.5, t = 7r/4, q = 3 and rout = 30. 


e is the emissivity of the material, rg is the radius of the 
emitting fluid element and fie is the cosine of the angle 
of the emission with respect to the disk normal. Eq. 
together with the delta-function expression for the line 
specific intensity gives 

dEo(Eo) = g'^€{re, fLe)S{Eo - gEint)dE ; (26) 

the extra factor of g in Eq. compared to Eq. comes 
from the change of argument in the delta-function. Now 
let Te and (jie be the plane polar coordinates in the disk 
of the emitting point (where = 0 is along the line 
of nodes where the disk intersects the observers plane 
of the sky). In addition let a and /? be the cartesian 
coordinates in the observer’s plane of the sky (where the 
X axis appears to lie along the lines (pe = Q and (j)e = tt), 
and To be the distance to the source. It follows that 
the solid angle element is given by d5 = dad/3/rg, and 
therefore integrating Eq. gives 

Fo(Eo) = ^ JJ dadfi g'^e{re, fie)5{EQ - gEi^^t) ■ (27) 

If light-bending is neglected then the polar coordinates of 
the emitting point in the disk (cg and (pe) may be related 
to the cartesian coordinates of the point in the observer’s 
plane of the sky (ag and /3g) via the transformation, 

Ue = Te cos (pe i Pe = Te Sin pe COS b , 
with Jacobian = r-g cos i , (28) 


where l is the inclination from which the disk is viewed 
(i.e. the angle between the observer’s line of sight and 
the spin axis of the BH). For further discussion of this as¬ 
sumption see Sec. |III C| With this in hand we may change 
the integral in Eq. |27| from being over the observer’s plane 
of the sky to being over the disk itself (Dropping the sub¬ 
script e), 


Fo{Eo) 



rdrdp g^e(r, fi)S (Eq - gE^nt) ■ 


(29) 

Hereafter the emissivity is taken to be a function of r 
only (no fie dependence), and is parameterised as a power 
law, e = where q will be referred to as the emissivity 
index. We have now reduced the problem to finding the 
red-shift factor as a function of the position in the disk. 
The red-shift factor is defined by 


9 = 


Eo 

Ee 




I observer 




emitter 


(30) 


where p^ is the four-momentum of the photon linking 
the emitter to the observer, is the four-velocity of the 
observer, and is the orbital four-velocity of the emitter 
(see Eq. 21). Since the metric is independent of both the 
t and p coordinates, the corresponding components of 
the photon’s four-momentum are conserved; 


[Pfj] = {Pt,Pr,Pe,P<j>) = {-E,pe,pe,A) ■ (31) 
















7 


If the observer is at infinity and is at rest with re¬ 
spect to the black hole the numerator of Eq. 32 is 


simply given by observer = —E. Using the orbital 
four-velocity from Eq. the denominator of Eq. is 
I emitter = —Eu^ + Qu*A, and hence the red-shift fac¬ 
tor is given by 


9 = 


d(l-UA) ’ 


(32) 


where A = A/E. Since A is conserved along the photon’s 
worldline it may be evaluated at infinity where the indices 
on the photon four-momentum may be raised and lowered 
using the usual Minkowski metric. At large distance the 
ratio of the contravariant components of the photon’s 
four-momentum gives the azimuthal impact parameter 


a, defined in Eq. 28 


a = -- 


rp^ 


= ■ ( 33 ) 


sm i 



the orbiting material is moving at relativistic speeds and 
for large inclinations material on one side of the disk is 
moving towards the observer while on the other side it 
is moving away. This produces both a blue and a red- 
shifted horn to the profile, the blue-shifted horn is always 
more intense due to the relativistic beaming effect which 
enhances the intensity of light emitted in the direction of 
travel (the headlight effect). 

The profile is also strongly sensitive to the radius of 
the inner edge of the disk, which in the Kerr case (under 
the modelling assumptions made here) is monotonically 
related to the spin parameter, a (see Eq. 181. For higher 
spins (and smaller values of risco) it can be seen from Fig. 
[^that the red-shifted wing of the line profile is much more 
prominent. Under the assumptions that the black hole 
is Kerr and the inner edge of the disk is located at the 
ISCO the line profile can be used as an accurate probe of 
the black hole spin, indeed such observations are among 
the best evidence for the existence of near maximally 
spinning black holes in some active galactic nuclei (see, 
for example, [20] )• 


The line profile also depends strongly on q; this pa¬ 
rameter must be fit for simultaneously with all other pa¬ 
rameters. The fiducial value q = 3 was chosen as this 
is approximately the form of the gravitational energy re¬ 
lease per unit area of the disk (see Fig. [^. Likewise the 
value of Tout must also be fit for, but this is less prob¬ 
lematic because (at least for rout 20) the profile only 
depends weakly Tout- 


B. Black-body spectra 


FIG. 2. The Flux as a function of radius, for various values of 
the spin parameter. Also shown is the power law, {r/M)~^. 
All other parameters were set to their fiducial values, l = 7 r/ 4 , 
M = Mq, Mo = lO“^^M 0 /yr and rout = 30. 


If light-bending is neglected then the photon linking 
the emitting point in the disk to the observer at infinity 
travels on a straight line with impact parameter a = 
rcoscj) (see Eq. 28). Substituting these results into Eq. 


32 gives the red-shift; 


/ j\ \/ iott + + Q pip^) 

g (r, (p) = -----—-. (34) 

I -h Ur cos 4> sm i 


The flux integral in Eq. [^may now be simply evaluated 
numerically. This was done initially for the unperturbed 
Kerr metric in order to reproduce known results (see, 
for example, Ea EHI) and examine the dependence of 
the spectra on various parameters. The results of these 
calculations are shown in Fig. For all calculations in 
this paper, unless otherwise stated, the disk parameters 
were set to the following fiducial values; a — 0.5, t = 7r/4, 
q = 3 and rout = 30. 

It can be seen from Fig. that the profile is strongly 
sensitive to disk inclination. Qualitatively this is because 


The line emission profiles calculated above encode how 
the gravitational and Doppler shifting (and lightbending 
if included) affect a single frequency source in the disk. 
If instead there is a broadband source the resulting spec¬ 
trum may be obtained by convolving the spectrum in the 
source’s rest frame with the line profile obtained in the 
preceding section. This method can be used to calculate 
the black body spectra by convolving with the Plank dis¬ 
tribution of a black body; this approach is described by 

m- 

Here we take the alternative, and more physically mo¬ 
tivated approach of directly integrating the flux over the 
disk to find the spectra. From Eq. [^ the radial flux 
(shown in Fig. is known, and assuming that the radi¬ 
ation is that of a black body the Stefan-Boltzmann law 
gives the radial temperature distribution in the disk 


r(r) 



(35) 


where a is the Stefan-Boltzmann constant. Hence the 
thermal spectrum may be written as the following inte¬ 
gral where the gravitational effects on the spectrum are 
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FIG. 3. The dependence of the Kerr black body spectrum on various parameters. Unless indicated otherwise all parameters 
were set to their hducial values, a = 0.5, b = 7 r/4, M — Mq, Mo ~ 10 and rout = 30. 


included through the redshift factor. 


Po{Eo) — 


8 cost 


TT 


r2Tr 


rd(j) dr 


^0 




• (36) 


As was the case for the line emission in Eq. the ther¬ 


mal spectrum in Eq. 36 is now in the form of an integral 
over the disk in r and (j) coordinates; this was evaluated 
numerically. Several example thermal spectra are shown 
for the Kerr metric in Eig. 

The thermal spectra of the disk is the convolution be¬ 
tween the line emission spectra and the Planck distribu¬ 
tion. Since the Planck distribution contains no informa¬ 
tion about the gravitational field the thermal spectra in 
Pig. [3] contain the same information about the black hole 
as the line spectra in Eig. only significantly smoothed 
out. Because of this smoothing it appears that the ther¬ 
mal spectra are less distinctive, and that therefore it 
would be harder to measure the parameters of the black 
hole using thermal emission than line emission. How¬ 
ever, the thermal spectra may be observed across a much 
wider frequency range and at a larger signal-to-noise ra¬ 
tio. Therefore it is not obvious a priori which method 
offers the best opportunity to constrain the metrics in 
Sec. |TT] In practice the most suitable technique depends 
upon the mass of the black hole; Iron line emission is typ¬ 
ically most suitable for supermassive black holes whilst 
thermal emission is used for the hotter disks around stel¬ 
lar mass black holes. In rare instances both techniques 
may be used simultaneosuly, and they have been demon¬ 
strated to give consistent resuts [4^ . Henceforth we con¬ 
sider only the Iron line emission, and we leave a detailed 
study of parameter estimation using thermal emission to 


future work. 


C. The effect of lightbending 


The formalism for calculating both the Iron line and 
thermal emission outlined in the preceding sections as¬ 
sumed that points in the image plane could be related 
to points in the disk by straight lines, in the sense that 
the Boyer-Lindquist-like coordinates were treated as if 
they were spherical polar coordinates in flat space (these 
assumptions are summarised in Eq. 28). As the light 
originates from the strong gravitational field the effects 
of lightbending (including frame dragging, if a 0) may 
be significant. 

On the other hand it may be hoped that the effect of 
lightbending will vary slowly with changing system pa¬ 


rameters, 9. If this is the case then while the inclusion 
of lightbending will have a significant impact on the ob¬ 
served profile it will have only a limited effect on our 
ability to measure the disk parameters, and hence on our 
ability to constrain deviations from the Kerr metric. In 
the language of the Fisher matrix calculations in Sec. |IV[ 
even if the spectra depends strongly on whether or not 
lightbending is included if the derivatives of the spec¬ 
tra do not then the fisher matrix remains unchanged. In 
order to test the validity of this assumption a small num¬ 
ber of calculations were performed including the effect of 
lightbending, and the results compared. 

A method for ray-tracing photons from the image plane 
to the disk was outlined in [41]. The method uses the 
fact that the spacetime is stationary and axisymmetric 
to write the t and cj) geodesic equations for the photon 
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FIG. 4. A set of images showing the resolved appearance of an accretion disk around a Kerr black hole. The colour indicates 
the redshift of the light from that portion of the disk; the top row shows the disk with no lightbending and the bottom row 
shows the disk including lightbending. The inclination angle is varried between plot; from left to right it takes the values 
t = {tt/IO, 27r/10, Stt/IO, dyr/lO}. The spin parameter and outer radius of the disk were fixed to a = 0.5 and rout = 10. 


in first order form, but integrates the r and 6 geodesic 
equations in second order form. This is ideal for our 
present purpose as the method does not require the ex¬ 
istence of a Carter-like constant, which does not exist in 
all the metrics discussed in Sec. im If a fourth and final 
constant of motion does exist (as is the case for the Bn 
metrics) then evaluating it along the resulting trajectory 
provides a convenient check on the numerical accuracy of 
the integration. 

The effect of lightbending on the appearance of a spa¬ 
tially resolved disk is shown in Fig. for varying inclina¬ 
tion (the colour scale indicates redshift factor, g). Partic¬ 
ularly at high inclinations the inclusion of lightbending 
significantly alters the appearance of the disk; the effect is 
most pronounced for large values of spin and inclination 
where the light from the far-side, inner edge passes very 
close to the horizon. However, the disk is not spatially 
resolved by our telescope, instead we observed the inte¬ 
grated flux across the disk. This shows a much smaller 
difference. The effect of including lightbending is also 
more significant for higher values of the spin parameter, 
because the disk extends closer to the black hole where 
the gravitational effects are stronger. 

For example, for the fiducial black hole parameters the 
Fisher matrix estimates an error on the dimensionless 
spin parameter of Aa = 0.06 if light bending is neglected 
(see right hand panel of Fig. [^. If instead lightbending 
is included then the same Fisher matrix analysis yields 
an error estimate for the spin parameter of Aa = 0.03. 
Changes in the error estimates of a factor of ~ 2 were 
observed for the other parameters. 


As anticipated above larger differences were obtained 
for higher values of spin and lower values of inclination, 
and smaller differences in the opposite limits. These 
changes are not enough to affect our conclusions in this 
paper. 


IV. THE FISHER MATRIX 


For each black hole metric described in Sec.|TT]the Iron 
line and thermal spectra are characterised by a small 
number of parameters 0; e.g. for the Iron line spectra 
h{0), 9 = (A,a,L,q,rout,^, where A is an overall multi¬ 
plicative factor relating to the unknown total luminosity 
and distance to the source, a is the spin parameter, l is 
the disk inclination, q is the emissivity index described 
in Sec. Ill A[ Tout is the outer radius of the disk and e is 
the vector of any metric deformation parameters. Given 
a measured spectrum, s, the challenge is to infer the pos¬ 
terior probability density on these parameters, P{9\s,I). 
The peak of this distribution is positioned at the best 
estimate of 9, and the characteristic width of the peak in 
each parameter direction indicates the uncertainty. The 
posterior probability density is related to the likelihood 
of the data given the parameters, P{s\9,I), via Bayes 
theorem. 


pm 


P{s\9)P{9) 

Pis) ’ 


(37) 
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where P{0) is the prior probability density of the param¬ 
eters and P{s) is the normalisation constant known as 
the evidence, 


P{s) = j AeP{s\e)P[9) . 


(38) 


For all calculations performed in this paper flat priors 
were assumed over all physically allowed regions of pa¬ 
rameter space so the posterior is simply proportional to 
the likelihood within this region. 

In order to calculate the likelihood it is first neces¬ 
sary to make some assumptions about the performance 
of our detector. For simplicity we assume that across the 
entire energy (frequency) range the detector has a con¬ 
stant energy (frequency) resolution /S.E (A/), and that 
in each energy (frequency) bin there is an independent 
Gaussian error of size a. The values of Ai? = 100 eV 
(A/ = 2.4 X 10^"^ Hz) and cr = 0.05 were choosen to 
make the hypothetical instrument broadly equivalent to 
the best current observations. For instance, the error on 
the spin measurement for the fiducial disk parameters 
under these assumption is Aa = 0.06 (see, for example, 
the results in Fig. [^. 

These simplifying assumptions neglect some important 
effects. For example, the frequency resolution of most 
X-ray detectors changes substantially across observable 
bandwidth, the random errors in each frequency bin are 
due to variations in the arrival rate of photons which is 
a Poisson not a Gaussian process (although in the limit 
of high signal-to-noise the assumption of Gaussian errors 
becomes correct), and in addition to the uncorrelated 
random errors there will be systematic errors which may 
be correlated between frequency bins. Nevertheless, un¬ 
der these simplifying assumptions the likelihood is given 
by the following expression. 


m= 


exp ( ^ ^ 


y/^J 


N 


exp l^s — h{9)\s — h{9)^'^ 


where the inner product has been defined as 


(39) 


{a\b) = 


^ u 

E O-i Oi 
i—1 


(40) 


In general the likelihood (Eq. pl9[ ) is a complicated func¬ 
tion of the parameters, 9. However, expanding the signal 
about the true parameter values, 9q, (using the Einstein 
summation convention) 


h{9) = h(0;) + ^ 


9=00 


69, + 0{{S9,f 


(41) 


and using s = n + h{9o) (where n is the particular reali¬ 


sation of the noise observed in the detector) gives 


c{9y 


exp (n|n) - 2 + P‘rj59i59j ^ 


where = 


/ dh{0) 




N 


9=00 


dh{9) 


89, 


0=00 , 


(42) 


Therefore, within the linear signal approximation used 
in Eq. the likelihood is a multivariate Gaussian, 
peaked a noise-realisation-dependent distance away from 
the true parameters, and with a covariance matrix given 
by the inverse of the Eisher information matrix, E^. An 
estimate for the error in each parameter may be read 
off from the corresponding component of the covariance 
matrix, 


A9,= 



(no sum on i). 


(43) 


The Fisher matrix formalism was used first to estimate 
the parameter estimation accuracy for a disk in the Kerr 
metric (Sec. VIA below) and subsequently for all of the 
bump y black hole spacetimes discussed in Sec. [H] (Secs. 
VI B| to VID below). In the case of the bumpy black 


hole spacetimes the true value of the bump parameters 
were set to zero, i.e. the Kerr metric was used, and the 
error estimate obtained for the bump is reported. The 
error on the bump parameter(s) are then interpreted as 
an estimate of the bound it may be possible to place on 
the size of the bump; i.e. if the true value of the bump 
parameter took this value then it would be marginally de¬ 
tectable with these observations. This bound should be 
interpreted as a lower limit, i.e. a best case scenario; in 
reality even if a non-zero value of a given deformation pa¬ 
rameter was returned in a particular experiment it would 
still be a non-trivial task to rule out more mundane ex¬ 
planations. For example, before claiming a detection of 
a deviation from the Kerr solution it would presumably 
be necessary to consider more complicated forms for the 
radial emissivity, e(r), than a simple power law. The free 
parameters in this new emissivity law would then have 
to be marginalised over and this would have the effect 
of increasing the errors on the other parameters. Other 
possibilities must also be considered; for example thick 
accretion disks, emitting material within the ISCO, re¬ 
processing of the light by surrounding material, etc. 

The applicability of the Fisher matrix rests on the va¬ 
lidity of the linear signal approximation in Eq. 41 this 
must hold at least within a few standard deviations from 
the peak in all directions in parameter space. In general 
it is impossible to know from the Eisher matrix alone 
whether one is within the region where the linear signal 
approximation may be safely applied. This question of 
the applicability of the Eisher matrix is addressed in Sec. 

IS 
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FIG. 5. A comparison of the likelihood surface predicted by the Fisher matrix (indicated by Icr (black) and 2a (red) contours) 
and the true posterior (shown as a density plot using samples produced from an MCMC algorithm). All parameters were set 
to their fiducial values; a = 0.5, l = 7r/4, q = 3 and rout = 30. 


V. VERIFICATION OF THE APPLICABILITY 
OF THE FISHER MATRIX FORMALISM 


The gold standard for parameter estimation is to nu¬ 
merically calculate the likelihood (or in general the pos¬ 
terior) surface over the region of parameter space of in¬ 
terest. This may be achieved in low dimensional prob¬ 
lems by simply evaluating the likelihood function on a 
grid of parameter points. Alternatively, and more effi¬ 
ciently in high dimensional problems, there exist a va¬ 
riety of Markov chain Monte Carlo (MCMC) algorithms 
designed to sample points from the target probability dis¬ 
tribution; the simplest of these is the Metropolis-Hastings 
algorithm. In Secs. |V A| a MCMC analysis is performed 
on a sample Iron line spectra for a typical case; the result¬ 
ing likelihood surfaces are compared with those predicted 
using the Fisher matrix. 


MCMC type algorithms quickly become (prohibitivly) 
expensive as the dimension of the parameter space in¬ 
creases. A consistancy check on the Fisher matrix was 
proposed by [42] which is relativly quick to implement. 
The consistancy check involves testing the validity of the 
linear signal approximation. For some central parame¬ 
ter values, 9q, a point 9 is picked at random from the la 
surface estimated by the Fisher matrix. The likelihood is 
then evaluated exactly, and approximated with the linear 
signal approximation, the ratio of the two likelihoods is 


denoted r[9). The logarithm of this ratio is given by 

logr(dl| = i - Ah(0l| - Ah(0l) , 

where A0 = 9 — 9q , Ah{9) = h{9) — h{9o ), 
dh{9) 


and hi = 


d9i 


(44) 


Small values of |logr(0)| incate that the linear signal 
approximation is holding out as far as the Icr surface 
in that particular parameter direction. This procedure 
may then be repeated for many points drawn randomly 
from the Icr surface to assess whether the approximation 
holds in all directions. It should be stressed that this 
only checks the internal consistency of the linear signal 
approximation and does not guarantee the accuracy of 
the Fisher matrix. In Sec. |VB| this consistency check was 
performed for the Iron line emission likelihood surface, as 
the consistency check is faster than a full MCMC it was 
performed for a range of spin and inclination parameter 
values. 

These two checks on the applicability of the Fisher 
matrix give us increased confidence in the results for 
the bounds on the various metric deformations found in 

Sec. ED 


A. MCMC 


A simple Metropolis-Hastings MCMC was used to 


sample from the likelihood distribution in Eq. 39 with 
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FIG. 6. Each panel of this figure gives the error in the column parameter versus the value of the row parameter for fixed fiducial 
values of the other parameters; a = 1, l = 7r/4, q = 3 and rout = 30. 


the Iron line spectra described in Sec. |III A[ The result¬ 
ing chain, plotted as a density histogram, is shown in 
Fig. [5] alongside the Gaussian contours from the Fisher 
matrix analysis. From Fig. [^it can be seen that there 
is unquestionably additional structure in the true poste¬ 
rior which is not captured by the Fisher matrix; however 
the widths of the Fisher matrix Gaussians give a good 
indiction of the scale of the true posterior. 


B. Mismatch ratio 


_ _7r_ _ 27r _ 371 _ 4?! 

^ ~ In^ ~ Tn^ ~ Tn^ ~ Tn 


a = 0.1 
a = 0.3 
a = 0.5 
a = 0.7 
a = 0.9 


0.0044/0.013 0.0088/0.032 0.025/0.15 0.10/0.50 

0.0046/0.016 0.0091/0.030 0.022/0.064 0.032/0.14 
0.0037/0.0069 0.011/0.026 0.019/0.083 0.028/0.077 
0.0039/0.0095 0.012/0.036 0.017/0.061 0.028/0.088 
0.0078/0.015 0.014/0.028 0.027/0.053 0.031/0.10 


TABLE I. All entries are of the form fJ,/M, where is the 
mean value of the logarithm of mismatch ratio calculated on 
100 points selected uniformly from the la surface and M is 
the maximum value of the mismatch ratio on the same set of 
points. All other parameters were set to their fiducial values 
of g = 3 and Tout = 30. (The forward slash between entries 
does not indicate division.) 


For a range of values of a and l the Fisher matrix was 
evaluated, and used to choose 100 points, randomly dis¬ 
tributed, on the Icr surface. The mismatch ratio was 
evaluated for all of these points and the results are sum¬ 
marised in Tab.[l] all values are less that unity, indicating 
that the Fisher matrix performs well for this problem. 


VI. RESULTS 

A. Kerr errors: Iron line 

Before attempting calculations in the bumpy black hole 
spacetimes the Fisher matrix formalism was used to as- 
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FIG. 7. The left-hand panel shows a pair of Iron line spectra for accretion disks around a Kerr black hole and a quadratic in 
spin CS black hole (CS2 metric) with a large deformation, C = 1- Also shown in the right-hand panel is the residual plot. All 
parameters were set to their fiducial values; a = 0.5, 4 = -tt/I, g = 3 and rout = 30. 


sess the accuracy with which it is possible to measure the 
standard parameters {A, a, t, q, X} with the Iron line ob¬ 
servations described in Sec. IIII Al The results are shown 
in Fig. In particular we note that the spin parameter 
can usually be measured with an error of Aa ~ 0.1. 

B. KS metric 

The formalism described in Sec. m was used to cal¬ 
culate the Iron line profile in the KS metric, Eq. [Tj As 
described in Sec. |III| the shift in the position of the ISCO, 
relative to the Kerr values, to first order in the small pa¬ 
rameter l/u! may be calculated; risco = + ^'^isco, 

where 

A -11 59o 

" 36u; 54^6.; ' ^ ^ 

The ISCO moves inwards for increasing deformation, 
therefore the effect of a large value of l/ui is to boost 
the redshifted red wing of the line profile relative to the 
blue-shifted peak, similar to the effect of increasing spin; 
see Fig. 

The Icr and 2 (t contours from a Fisher matrix analysis 
are plotted in Fig. it can be seen that there is a rather 
stark degeneracy between the spin and deformation pa¬ 
rameter. By comparing Fig. |9]to Fig.j^it can be seen that 
the errors in all other parameters are virtually unaffected 
by the inclusion of the deformation parameter. The effect 
of the degeneracy, apart from inhibiting any measure¬ 
ment of the spin parameter, is to make it very difficult 
to place any bound on the deformation; Tab. [TT| gives the 
bound it is possible to place on oj for different values of 
spin. It should be remembered that the KS metric is 
valid only to linear order in a, and in addition the metric 
exhibits unphysical properties (e.g. naked singularities, 
closed timelike curves, etc) for a; < 1/2. Therefore, bear¬ 
ing in mind the fact that the bounds derived from the 
method discussed in Sec. should be treated as lower 


Redshiftz 

1 0.33 0 - 0.2 - 0.33 



Frequency/ fo 

FIG. 8. A series of Iron line spectra for accretion disks around 
KS black holes varying the value of to. All other parameters 
were set to fiducial values, a = 0.5, 4 = 7r/4, g = 3 and 
rout = 30. 


limits, the conclusion to be drawn from Tab. [IT] is that it 
would be extremely difficult to place any constraint on 
the KS deformation parameter using observations of this 
type. 

There have been several attempts to place bounds on 
the KS deformation parameter using solar system tests of 
gravity such as perihelion precession, deflection of light 
by the Sun and radar echo delay observations 
The bounds obtained from these studies are typically on 
the length scale w > 5 x 10“^® m, corresponding to a di¬ 
mensionless bound ^ 3 x 10“^®, with M = Mq. It 
should be noted that these bounds are much less strin¬ 
gent than the requirement imposed here that uj > 1/2, 
which is necessary to ensure that there is an event hori¬ 
zon. However, as these tests were conducted in the weak 
field around a material object whose radius is much larger 
than the gravitational radius the vacuum solution is not 
valid down near the event horizon and the constriant 
w > 1/2 need not apply. 
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FIG. 9. The la and 2a contours from a Fisher matrix analysis on the KS metric. It can be seen that introducing the extra 
degree of freedom in the KS deformation parameter has introduced a clear degeneracy with the spin parameter. From the 
bottom righthand plot it can be seen that the data is almost equally consistent with any value of cu in the range (1/2, oo), so 
no bound may be placed in the deformation. All parameters were set to their fiducial values; a = 0.5, t = vr/l, 5 = 3, rout = 30 
and Y = l/w = 0. 


C. CS metric 


Shown in Fig. is the Iron line profile for a Kerr black 
hole and a CS2 black hole with deformation C = 1. It 
should be remembered that the CSl and CS2 metrics 
in Eqs. and IIA 3 are only valid in the limit C ‘C 1. 
However, even with the large deformation used in Fig.[^ 
the change in the line profile is only visible in the residual 
plot shown in the right-hand panel. The Iron line profile 
for the CSl metric shows very similar behaviour to that 
in Fig.0 


Bearing in mind the results in Fig. it is clear that 
it will be extremely difficult to bound the CSl deforma¬ 
tion parameter using this technique. Tab. |III| shows the 


bounds it is possible to place on both the CSl and CS2 
deformation parameters for a range of values of a; it was 
found that no bounds less than unity were possible with 
Iron line observations, however, the best results were ob¬ 
tained for the CSl metric and high values of spin. Both 
the CSl and CS2 metrics are expansions in the a param¬ 
eter, therefore the bounds for the higher values of spin 
(particularly a = 0.9) should be treated with more cau¬ 
tion than the low spins. However, this does not effect 
our main conclusion that no bounds less than unity were 
possible. 

For comparison, weak field tests using the frame¬ 
dragging effect around the Earth measured by the Grav¬ 
ity Probe B and the LAGEOS satellites places a bound 
< 10® km [30]. Tighter bounds will be possible using 
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A(l/tj) 

a = 0.1 

15.2 

a = 0.3 

9.23 

a = 0.5 

5.47 

a = 0.7 

4.76 

a = 0.9 

3.44 


TABLE II. The Bounds it is possible to place on the small KS 
deformation parameter (l/w) for different values of the spin 
parameter, a. All other parameters were set to their fiducial 
values; b = tt/A, q = 3 and rout = 30 and Y = 1/cu = 0. 


strong field tests, for example it was found that bounds 
of km corresponding to a dimensionless bound 

of C < X10“^ would be possible with eLISA observations 
of EMRIs [47]. 



ACcsi ACcS2 

a = 0.1 

101 122 

a = 0.3 

97.0 62.6 

a = 0.5 

24.0 25.4 

a = 0.7 

23.9 9.98 

a = 0.9 

1.32 5.34 


TABLE III. The Bounds it is possible to place on the CSl 
and CS2 deformation parameters for different values of the 
spin parameter, a. It should be remembered that the CS 
parameter is constrained to be ^ <C 1, therefore no meaningful 
constraint may be placed with these observations. All other 
parameters were set to fiducial values, b = it/ A, q = 3, rout = 
30 and C = 0. 


D. Bn metrics 


Shown in Fig. 10 are a series of Iron line profiles for the 
metrics defined in Sec. jlT] with the constants B 2 varied. 
The deformation parameters were varied between 0 and 
0.5. Shown in Fig. [^is the bounds it is possible to place 
on the different deformations using Iron line observations 
for different values of a. 


The results in Fig.fTTjshow that the tightest constraints 
on the Bn bumpy black holes can be placed on the low¬ 
est values of N for the highest values of spin. The lower 
values of N have deformations entering at lower powers 
of 1/r, since all the emission originates from outside the 
horizon (where r > 1) it is to be expected that defor¬ 
mations at lower N are easier to constrain. It is also to 
be expected that higher values of spin make placing con¬ 
straints easier, because the ISCO moves to smaller values 
of r for larger spins, and most of the emission comes from 
close to the ISCO, the deformation has a greater effect 
on the spectra for high spins. 


VII. DISCUSSION 

Observations of accretion disks offer an enticing oppor¬ 
tunity to probe strong gravitational fields. Such observa¬ 
tions have already provided some of the best evidence for 
the existence of black holes (in particular the existance 
of an event horizon, and a separate innermost stable cir¬ 
cular orbit) and are routinely used to measure the spin of 
black holes. In this paper we have considered the possi¬ 
bility of using X-ray emission from a thin accretion disk 
to distinguish between the Kerr black hole predicted by 
GR and alternative black hole metrics. For this purpose a 
series of parametrically deformed Kerr metrics, referred 
to as bumpy black holes, expanded in powers of M/r, 
were used [221 [28] . In this paper both the transition line 
and the thermal emission around these general bumpy 
black holes were calculated. The Fisher matrix formal¬ 
ism was used for the Iron line emission to estimate the 
accuracy to which the various disk parameters could be 
measured. 

In addition to this general family of bumpy black holes 
the disk emission in several specific, known black hole 
spacetimes were also considered. In particular the Iron 
line emission in the metric CS2 due to [23] was considered 
for the first time in the literature. For the CS2 metric it 
was found that it is impossible to place constraints on the 
small, dimensionless coupling parameter less than unity. 

The technique of using accretion disk observations to 
constrain alternative black hole solutions has limitations. 
Qualitatively, the method works well for measuring the 
spin of a Kerr black hole, but adding extra deformation 
parameters typically introduces degeneracies which make 
constraining the bump parameters very difficult. This 
was very clearly seen in the case of the Kehagias Sfetsos 
metric, see Fig. [^ For the BN metrics it was found that 
bounds less than unity were only possible in a limited 
number of cases, typically for the leading order deforma¬ 
tions around highly spinning black holes. It should also 
be remembered that the bounds presented here are opti¬ 
mistic, “best case” scenarios for the reasons discussed in 

Sec.ia 

In the near future it will become possible to perform 
this sort of test of GR using gravitational wave obser¬ 
vations of the inspiral and coalescence of a binary con¬ 
taining neutron stars or black holes. In particular, the 
possibility of using eLISA observations of extreme mass 
ratio inspirals to constrain the GSl metric was recently 
considered in m- It was found that bounds on the di¬ 
mensionless deformation parameter e < 10“® may be pos¬ 
sible, several orders of magnitude better than the bounds 
estimated here. As the GSl metric is a special case of the 
general Bn metrics considered here, it would be interest¬ 
ing to try to constrain the Bn metrics here using the same 
techniques and see if similar improvements are possible. 
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FIG. 10. A series of Iron line spectra for accretion disks aronnd B2 bumpy black holes with varying values of the bump 
parameters. All other parameters were set to fiducial values; a = 0.5, t = tt/4, g = 3, rout = 30, and yij = 0 unless otherwise 
indicated. 
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FIG. 11. The bounds it is possible to place on the different Bn bnmp parameters for different valnes of the spin parameter, 
a. The tightest bonnds can typically be placed for low valnes of N (i.e. B 2 ) and for the most highly spinning black holes. All 
other parameters were set to fiducial values; i = n/4, q = 3, rout = 30 and e = 0. 
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Appendix A: Metric components 


2. CS2 


1. BN 

For the bumpy Bm metrics discussed in Sec.|n]all of the 
metric coefficients up to 0{\/r^) are reproduced here. 

htt ,2 = 7i,2 + 274,2 - 2073,1 sin^ 6 
htt ,3 = 71,3 - 874,2 - 271,2 + 274,3 + 8073,1 sin^ 9 
huA = -874,3 - 271,3 + 274,4 + 874,2 + 7i,4 - 8073,1 sin^ 9 
(71,2 + 274,2) sin^ 9 + 20^73,1 cos^ 9 sin^ 9 
htt,5 = 160^73,1 sin'^ 9 + sin^ 9 [4073,3 + o^ (71,3 - 271,2 
— 1274,2 + 274,3) — 12o^73,i] + o^ (874,2 + 271,2) 
+71,5 + 274,5 — 271,4 + 874,3 — 874,4 (Al) 


ht 4,,2 = -M sin^ 9 [ 73,3 + o ( 71,2 + 74 , 2 ) + a^ 73 .i] 

ht<p,3 = sin"* ^ 

+M sin^ 9 [(273,3 - 73,4) + a (674,2 - 74,3 

+ 271,2 — 71 , 3 ) + 273 , 10 ^] 

ht<t,A = M sin"* 9 [o^ (873,1 - 73,3) + (-71.3 - 74,2) 

-0^73,1] + sin^ 9 [(273,4 - 73,5) + o (-74,4 
-874,2 + 674,3 - 71,4 + 271,3) - 0^73,3] 
ht,p,5 = -16Mo'‘73,i sin® 9 

+M sin"^ 9 [a^ (-273,3 - 73,4) 

+a® ( 74,3 + 1074,2 + 271,2 - 71 , 3 ) + 140^^73, 1 ] 

+ sin^ 6> [(273,5 - 73 , 6 ) 

+a (-71,5 - 874,3 - 74,5 + 271,4 + 674,4) 
-73,40^ + o® (-271,2 - 674,2) - 2 o^73,i] (A 2 ) 


hrr ,2 — —71.2 

hrr ,3 = —71,3 — 271,2, 

hrrA = -71,4 - 271,3 - 471,2 + (1/2)71,20^ (1 - cos 20) 
hrr ,5 = sin^ 9 ( 71,3 + 271 , 2 ) 

-71,5 - 271,4 - 471,3 - 871,2 + 2o^7i,2 (A3) 


For the CS2 metric discussed in Sec. Ill A 3l the metric 
perturbations are reproduced here, with /(r) = 1 — ( 2 /r). 
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= 0 
! — 0 

^ 00,0 = 2M^a73,i sin'^ 0 
— 0 

sin'^ 0 [ 2073,3 + 0 ^ 71,2 + o® 73 .i (4-2 cos^ 0)] 
^ 00,3 = sin® 0 + sin'‘ 0 [a (- 473,3 + 273 , 4 ) 

+a^ (- 271,2 - 474,2 + 71 , 3 ) - 4 a® 73 ,i] (A4) 
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Appendix B: Radial dependence of the flux 


[33] derive various expressions for the radial structure 
of the disk, including an expression for the radial depen¬ 
dence of the flux (Eq. in the main text). Here the 
derivation is summarised for completeness. 

The thin disk is assumed to be axisymmetric, station¬ 
ary, and lying in the equatorial place; therefore all quan¬ 
tities in the disk depend only on the radial coordinate. 
It is assumed that the material in the disk moves (very 
nearly) on circular geodesics. The four-velocity of indi¬ 
vidual fluid elements, Uq, when mass averaged over the 
disk structure must therefore be the four-velocity of the 


geodesic orbit in Eq. 22 


1 


r-\-h p-\-h 

/ dz PqUq , where S(r) = / dz 

J-h J-h 


nr) J-H 

(Bl) 

Where po rnass density in the rest frame of the orbiting 
material. Without loss of generality the stress-energy 
tensor may be decomposed by writing 


T^^'' = po (1 + n) -f -t 2u^rq<') 


(B2) 


where the physical interpretation of each term becomes 
clear in the rest frame of the orbiting material; H is the 
specific internal energy, is the stress tensor in the 
averaged rest frame of the material and is the energy 
flow vector. The round brackets in the superscript of Eq. 
B2 denote symmeterisation with respect to the enclosed 
indices. The tensors and obey the orthogonality 
relations = 0 a nd = 0. By analogy 

with E(r) in Eq. Bl the average stresses in the disk are 
defined as 




w; = 


dzt' 


(B3) 


f —h 


much less than the acceleration due to gravity otherwise 
the material would be pushed off its geodesic trajectory. 
Using the approximate relation « poH (which is valid 
for any astrophysically reasonable matter |53|) leads di¬ 
rectly to the condition of negligible specific heat; 

radial pressure acceleration Ri |i9rUr/po| ~ \dr (Ur/po)| 
gravitational acceleration « \drE\ ~ \dr{l — E)\ 
\dr.{trr/po)\<^\dr{l- E)\ ^ H < 1 - E . (B6) 

If the internal energy is small compared to the gravita¬ 
tional potential energy, this means that as the material 
spirals in towards the black hole all of the gravitational 
potential energy is radiated away. 

With these simplifying assumptions in place the equa¬ 
tions governing the structure of the disk follow from the 
conservation of stress-energy = 0), and the con¬ 

servation of rest mass of the fluid [36] . 

{pour) = 0 . (B7) 


This is integrated over the spacetime volume 
{V : t e (to, to+T), r e (r, r -I- Ar), € (0, 27r), 

z G {—h, +h)}. Gauss’s theorem is then used to convert 
the volume integral into a surface integral over the 
boundary dV with area element Id^Al. 


0 = 


0 = Idv Pourn^\d^A\ 

I-H drdMz V=SPou* 

lo^ I-h dtd(j)dz yf^pou^ 


(B8) 


■ t - tQ-\-T 

. t^to 
r'—r-\-Ar 


I^H dtdrdz V^pou^ 


r'—r 

-| ^> — 277 


dtdrdc^ n^pou^ 


■ (p—O 

- z—-\-h 

( 

- z——h 


(B9) 


It is also assumed that heat flow within the disk is 
negligible except for in the vertical direction, which is 
reasonable as the disk is thin. The non-local heating 
effects due to light emitted by one portion of the disk 
being re-absorbed by another portion are also neglected. 


q* = = 0 , at z = ±h. 


(B4) 


Since the only time-averaged stress that reaches out of 
the disk to infinity is carried by photons (neglecting grav¬ 
itational radiation and any coherent superposition of long 
wavelength radiation), and using Eq. B4 on the upper 
and lower edges of the disk the following terms of the 
stress-tensor disappear; 


n = n =n=d 

and \q^\ = F{r ), at z = ±h . 


(B5) 


The first and third terms in the above expression are 
zero by the assumed stationarity and axisymmetry of the 
system. The final term is also zero because there is no 
motion in the vertical direction, = 0. Therefore Eq. 
|B9| simplifies to 

0 = 2'KT/S.r (y^S(r)u'’) ,, 

^ Mq = —27ri/—gE(r)u’' = constant, (BIO) 

where Mq is the accretion rate. 

The second conservation law is that of angular momen¬ 
tum. Again the differential form of the conservation law 
is integrated over the volume V and Gauss’s law used to 
turn this into a surface integral over dV. 

0 = V.jr where jr = 

\d(l) 

0 = I \d^A\ 

JdV 



Since we are assuming the particles in the disk are, very 
nearly, on circular, geodesic orbits it follows that the ac¬ 
celeration due to pressure gradients in the disk must be 


(Bll) 
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0 = 




dM0dz [po(l + + q 4 >vJ'] 


r'—r-\-Ar 



z—-\-h 


[po(l + n)it0M^ + + u^q^ + qci,u^] 


z— — h 


(B12) 


The (/) index has been lowered and the t and </> integral 
terms have been set to zero due to the assumed stationar- 
ity and axisymmetry of the system. Using the negligible 
internal energy condition derived above, Eqs. |B4| and |B5[ 
and the fact that = 0 this becomes 


0 = 


p-\-h 

2ttT / dz [poU 4 ,u^ + t^]^/^ 
J-h 


r'—r-\-Ar 


+ 


nr-\-Ar 

2ttT / dr u^q^yf^ 


z—-\-h 


z——h 


(B13) 


An extra factor of two has appeared on the left-hand side 


of Eq. B14 from the fact that a flux F{r) is radiated from 


both sides of the disk. 

The third and final conservation law is that of conser¬ 
vation of energy. By performing the same type of ma¬ 
nipulations to this equation as was done for Eq. |B12| we 
obtain, 

0 = where S'" =(^)^ , (B15) 

MoE + 2 -k./=IwA = 47ry=iE(r)£; . (B16) 


Making use of the orthogonality = 0 which implies 
that = 0 ^ W/ + = 0, this equation may 

be rewritten in terms of 
angular momentum equation. 


as was the case with the 


MqE - 27ry=^lU/U = 47ry=iF(r)U (B17) 


From Eqs. and [T7] it can be seen that the en¬ 

ergy, angular momentum and angular velocity satisfy the 
following energy angular momentum relation. 


E r = flT . 


(B18) 


Eqs. B14 and |B17 may now be integrated to find the 
radial dependence of the flux. This is done by multiplying 
Eq. |B14| by U and subtracting the result from Eq. |B17| 
to obtain an expression for the torque; 


gives a differential equation for F(r), 


(E-VLlf 

- F{r) 


n, 


= Mo{e- L^r, 


(B20) 

which may be readily integrated. To fix the constant of 
integration we use the zero torque boundary condition at 
the inner edge of the disk, F(risco) = 0. Therefore we 
have an expression for the radial flux from the disk. 


47rV-gF(r)Z = [MqA - 27rV-glE/j . (B14) F{r) = 




47rA/=g (^E - 


(^E - 


L rdr. 


(B21) 

For completeness we also present the final radial struc¬ 
ture expression derived in [^, the expression for the 
torque per unit circumference as a function of radius, 

WL 


w;{r)= 


— E — 
27rv/=^(E-ffL)' 

/ T 

{e - VlL^ L^rdr . 


(B22) 


W; = 2E{r) 


ni- E 


(B19) 


Substituting this back into Eq. |B14| and rearranging and 
using the energy angular momentum relation in Eq. B18 
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